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The self consistency between the impressive DAMA annual modulation signal and the 
differential energy spectrum is an important test for dark matter candidates. Mirror 
matter-type dark matter passes this test while other dark matter candidates, including 
standard (spin-independent) WIMPs and mini-electric charged particle dark matter, 
do not do so well. We argue that the unique properties of mirror matter-type dark 
matter seem to be just those required to fully explain the data, suggesting that the 
dark matter problem has finally been solved. 
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Recently it was pointed out[l] that the impressive DAMA/Nal annual modulation 
signal[2, 3] can be simply explained if dark matter is identified with mirror matter. 
The purpose of this article is to explore further this explanation of the DAMA/Nal 
experiment. 

Recall, mirror matter is predicted to exist if nature exhibits an exact unbroken 
mirror symmetry [5] (for reviews and more complete set of references, see Ref.[6]). For 
each type of ordinary particle (electron, quark, photon etc) there is a mirror partner 
(mirror electron, mirror quark, mirror photon etc), of the same mass. The two sets 
of particles form parallel sectors each with gauge symmetry G (where G = SU{3) ® 
SU {2)®U{1) in the simplest case) so that the full gauge group is G ® G. The unbroken 
mirror symmetry maps x —x as well as ordinary particles into mirror particles. 
Exact unbroken time reversal symmetry also exists, with standard CPT identified as 
the product of exact T and exact P[5]. 

Mirror matter is a rather obvious candidate for the non-baryonic dark matter in 
the Universe because: 

• It is well motivated from fundamental physics since it is required to exist if parity 
and time reversal symmetries are exact, unbroken symmetries of nature. 

• It is necessarily dark and stable. Mirror baryons have the same lifetime as ordi- 
nary baryons and couple to mirror photons instead of ordinary photons. 

• Mirror matter can provide a suitable framework for which to understand the large 
scale structure of the Universe[7]. 

• Recent observations from WMAP[8] and other experiments suggest that the cos- 
mic abundance of non-baryonic dark matter is of the same order of magnitude as 
ordinary matter fldark ^ ^b- A result which can naturally occur if dark matter 
is identified with mirror matter [9]. 

Ordinary and mirror particles interact with each other by gravity and via the 
photon-mirror photon kinetic mixing interaction: 

^ = l^'^^^;. (1) 

where F^^'^ {F'nv) is the field strength tensor for electromagnetism (mirror electromag- 
netism) ^. Photon- mirror photon mixing causes mirror charged particles to couple to 
ordinary photons with a small effective electric charge, ee[5, 11, 12]. 

Interestingly, the existence of photon-mirror photon kinetic mixing allows mirror 
matter-type dark matter to explain a number of puzzling observations, including the 
Pioneer spacecraft anomaly [13, 14], anomalous meteorite events[15, 16] and the unex- 
pectedly low number of small craters on the asteroid 433 Eros[17, 18]. It turns out 
that these explanations and other experimental constraints[19, 20, 21] suggest that e 
is in the range 

10-9 - |e| ~ 5 X 10"^. (2) 

^Givcn the constraints of gauge invariance, renomalizability and mirror symmetry it turns out[5] 
that the only allowed non-gravitational interactions connecting the ordinary particles with the mirror 
particles are via photon-mirror photon kinetic mixing and via a Higgs-mirror Higgs quartic coupling, 
C = X<j)^ (j)' . If neutrinos have mass, then ordinary - mirror neutrino oscillations may also occur[10]. 
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However, the strongest direct experimental evidence for mirror matter-type dark matter 
comes from the impressive DAMA/Nal experiment [2, 3]. This experiment, located 1.5 
km underground in the Gran Sasso Laboratory, is designed to directly detect dark 
matter particles using a Nal target. This experiment has been in operation for about 
a decade and the DAMA collaboration have claimed discovery of dark matter. The 
approach of the DAMA/Nal experiment is to measure nuclear recoils of Na, I atoms 
due to the interactions of dark matter particles. Because of the Earth's motion around 
the sun, the interaction rate should experience a small annual modulation: 



According to the DAMA analysis[2], they indeed find such a modulation over 7 annual 
cycles at more than 6a C.L. (in the 2-6 keVee energy range): 



Furthermore, their data is beautifully self consistent, because their data fit gives 
T — 1.00 ± 0.01 years and to — 144 ± 22 days, consistent with the expected val- 
ues. [The expected value for to is 152 days (2 June), where the Earth's velocity reaches 
a maximum with respect to the galaxy]. This self consistency is highly non-trivial: 
there is simply no reason why their data should contain a periodic modulation or why 
it should peak near June 2. In fact, the estimated systematic errors are several orders 
of magnitude too small to account for the signal[3, 4]. It seems that DAMA have 
indeed discovered dark matter as they have claimed. 

It was pointed out in Ref.[l] that the DAMA experiment is sensitive to mirror 
matter-type dark matter. Halo mirror atoms can elastically scatter off ordinary atoms 
via the photon-mirror photon kinetic mixing interaction, Eq.(l). The DAMA exper- 
iment is not particularly sensitive to very light dark matter particles such as mirror 
hydrogen and mirror helium. Impacts of these atoms (typically) do not transfer enough 
energy to give a signal above the detection threshold[l]. The next most abundant el- 
ement is expected to be mirror oxygen (and nearby elements). A small mirror iron 
component is also possible. Interpreting the DAMA annual modulation signal in terms 
of O', Fe', we found that[l]: 



where the errors denote a 3 sigma allowed range and ^a' = Pa'/{0-^ GeV/cm^) is the 
A' proportion (by mass) of the halo dark matter. 

The value of e suggested by the DAMA experiment, Eq.(5), would also have im- 
plications for the orthopositronium systcm[19]. The current experimental situation, 
summarized in Ref.[20, 21], implies that |e| ~ 5 x 10"'', which is easily consistent with, 
Eq.(5). Importantly, a new orthopositronium experiment has been proposed[21] which 
can potentially cover much of the e parameter space suggested by the DAMA experi- 
ment. Such an experiment is very important - not just as a check of the mirror matter 
explanation - but also because dark matter experiments are sensitive to e^/^A' and an 
independent measurement of e would allow the extraction of ^a' values. 



Acos27i{t -to)/T. 



(3) 



A = 0.019 ± 0.003 cpd/kg/keV. 



(4) 




(5) 
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In Ref.[l] we found that a DAMA/Nal annual modulation signal dominated by 
the Fe' component, is experimentally disfavoured for three independent reasons: a) 
it predicts a differential energy spectrum rate larger than the measured DAMA/Nal 
rate, b) potentially leads to a significant diurnal effect c) should have been observed 
in the CDMS experiment. Thus it seems probable that the mirror oxygen component 
dominates the DAMA annual modulation signal, which from Eq.(5) means that ^o' ~ 
4,^Fe'- In this case there are no significant problems. 

If the DAMA signal is dominated by O', then things depend on only one parameter, 
e\/^0'- This parameter is fixed from the annual modulation signal, Eq.(5), giving a 
definite prediction for the differential energy spectrum in the DAMA experiment [up 
to the 15% 1-sigma experimental uncertainty in the annual modulation amplitude, 
A, Eq.(4)]. The prediction is shown in Figure 1, along with the DAMA/Nal data 
(obtained from figure 1 of Ref. [22])^. Note that for the annual modulation analysis, 
the DAMA collaboration have cautiously adopted a 2 kcV software energy threshold. 
Nevertheless, for the differential energy spectrum, they have published data going 
down to 0.8 keVee[22]. In general I support caution, however on this occasion we will 
be adventurous; we will use the published data without further reservations, with the 
aim of pushing existing information to the limit. 




1 2 3 4 5 6 7 8 

Energy [keVee] 

Figure 1: DAMA/Nal differential energy spectrum prediction for He' dominated halo with 
a small O' subcomponent (with vo{He') = 2vo{0') = ^ km/s). The value of e-y/fo^ is 
fixed by the anmial modulation signal, with A = 0.019 cpd/kg/keV (central value) solid line, 
A = 0.022 cpd/kg/kcV (1-sigma upper limit) upper dashed line, A = 0.016 cpd/kg/keV 
(1-sigma lower limit) bottom dashed line. Also shown is the DAMA/Nal data[22]. 

^See Ref.[l] for details of the cross section, form factors and rate equations used (including effects 
of energy resolution, quenching factors etc). 
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The fit of the prediction with the experiment is striking for low recoil energies, as shown 
in Figure 2. 




1.4 1.6 
Energy [keVee] 

Figure 2: Same as Figure 1, but low recoil energy region magnified. 



It suggests a background contribution which is roughly constant above 3 keV and fades 
to zero below 2 keV, which seems plausible given the dip in the data at ~ 2 keV. 

Our analysis has assumed that the halo has the standard isothermal Maxwellian 
distribution (see the appendix for details): 



(6) 



In the case of a halo with mass dominated by one component, which we have assumed 
is He' (which is suggested by mirror BBN arguments [7]), then 



voiHe') 



Vrot 



220 



V3 V3 



km/s. 



(7) 



Even if the mass of the halo is assumed to be dominated by He' , there are, of course, 
still important astrophysical uncertainties here. We consider a range for vo{He') of 
(c.f. Ref.[23]): 

170 km/s ~ VSvo{He) ~ 270 km/s. (8) 



*In the earlier paper, Ref. [1], we set vo{He') ~ Vrot ~ 220 km/s. However, this neglects the 
pressure due to the mirror electrons, which are important since He' should be fully ionized in the 
halo, given that T keV (see the appendix for details). 
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Since the mirror ions/electrons are in thermal equihbrium with each other, knowledge 
of Vq for He' will fix vq for the other elements: 



vo{A')=vo{He')JMHe'/M^ 



(9) 



We now study the effect of varying vq on the mirror matter DAMA prediction, first 

by retaining the assumption of He' dominant mass halo and considering the variation 
of Voi^He') given in Eq.(8). In this case, y/3vo{He') = 220±50 km/s implies V3vo{0') = 
110±25 km/s. In going from \/3vo{He') = 170 km/s to \/3vo{He') = 270 km/s we find 
that the central value for e^/^ from the annual modulation signal changes by only 
~ 10%. In Figure 3 we examine the effect of varying vo on the predicted differential 
energy spectrum (while keeping e^/^o' fixed by the annual modulation signal, A = 0.016 
cpd/kg/keV). As the figure shows, the effect is very slight, with negligible modification 
to the predicted rate. 
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Figure 3: DAMA/Nal differential energy spectrum prediction, assuming tliat tlie mass of 
the halo is dominated by He' with a small O' subcomponent. The value of eV^o' is 
fixed by the DAMA annual modulation signal (taking the 1-sigma lower value, A = 0.016 
cpd/kg/keV), with V3vo{He') = 220 km/s (solid line), V3vo{He') = 170 km/s (dashed line) 
and VSvo{He') = 270 km/s (dotted line). 



A more extreme variation in vo{0') is possible if we vary the dominant halo component. 
In particular, consider the case where the mass of the halo is dominated by H' rather 
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than He' . This would mean that \/2vq{H') ~ Vrot ~ 220 km/s (see appendix) and 

v^{0') = Vo{H')^Mh'/Mo' ^ km/s. (10) 

This case, along with the extreme case of a pure O' halo, with vq{0') = 220 km/s 
(which we give just for curiousity) is shown in Figure 4. 



14 -\ 
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Figure 4: Examing the effect of varying the dominate halo component. In each case we 
assume that the halo contains a mirror oxygen subcomponent but the mass of the halo is 
dominated by a) mirror helium (solid line) b) mirror hydrogen (dashed line) and c) containing 
only mirror oxygen (dotted line, top). In each case the value of €\/^o' is fixed by the annual 
modulation signal (A = 0.016 cpd/kg/keV). 

Although the data seems to prefer a dominant He' halo, the case of a dominant H' halo 
is certainly still possible, after taking into account possible uncertainties. However, a 
pure O' halo appears to be completely ruled out. It is of course, theoretically highly 
unlikely too, but it does illustrate that the nice fit of the He' case (and to some extent, 
the H' case) with small O' subcomponent is quite non-trivial. It relics on the effect 
of the dominant light elements {He' and/or H' and the mirror electrons) to make 
vo{0') < 220 km/s. 

So far we have examined mirror matter dark matter because it seems to have the 
right properties to explain the experimental data. It also seems to be theoretically 
favoured, since it arises from just one fundamental hypothesis: that is the particle in- 
teractions respect an exact, unbroken mirror symmetry. We now briefly examine alter- 
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native, but less elegant (from a particle physics point of view!) dark matter candidates. 
In particular imagine that the non-baryonic dark matter is composed of mini-charged 
particles (and antiparticlcs) of electric charge ±ee and mass Mx- This case may look 
superficially similar to the mirror matter case, but with more freedom since the mass, 
Mx, can be arbitrary. However, since in this case the halo would have only one mass 
component we would expect Vq{X) ^ 220 km/s. Again we can fit the DAMA annual 
modulation signal, but such a model is not consistent with the measured differential 
energy spectrum. We illustrate this in Figure 5. This clearly demonstrates the ne- 
cessity of a multi-component halo to explain the data - which naturally occurs in the 
mirror matter model. 
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Figure 5: DAMA/Nal differential energy spectrum prediction for mini-electric charged par- 
ticle dark matter (assuming standard halo with vq = 220 km/s), corresponding to the three 
cases a) Mx = 40 GeV (dotted line), b) Mx = 20 GeV (dashed line) and c) Mx = 12 
GeV (solid line). The latter value is actually the best fit to the measured differential energy 
spectrum [for 5 < Mx/GeV < 100]. 



Finally in Figure 6 we compare mirror matter dark matter with standard spin 
independent (SI) WIMP dark matter. Recall WIMP dark matter is quite different 
from mirror matter dark matter because its cross section is point-like: 

da 27re2a2^2^'2 

Fj^[Er)Fj^i[Er) Mnror matter case 



da ^ 2G|MaAV^2 
dEn 7rv2 



(Er) wimp case (11) 
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where Ma is the target nucleus mass (with mass number ^4), g is the dimensionless 
couphng constant characterising the WIMP-nucleus interaction, v is the impact velocity 
(in lab frame) and the F's are the form factors. 
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Figure 6: DAMA/Nal differential energy spectrum prediction assuming standard (SI) WIMP 
dark matter. Three representative cases are considered a) 30 GeV WIMP mass (solid line) b) 
50 GeV WIMP mass (dashed line) and c) 100 GcV WIMP mass (dotted line). In each case 
the cross section is fixed to the value given by the best fit of the annual modulation signal 
{A = 0.019 cpd/kg/keV). We have assumed a standard halo with vo = 220 km/s. 



From Figure 6 wc sec that after fixing the parameters to the DAMA annual modu- 
lation signal the (SI) WlMPs arc only marginally consistent with the measured differ- 
ential energy spectrum and only if they are hght^. Note though that we have assumed 
a standard halo with vq = 220 km/s. Lower values of Vq and non-standard halo models 
can improve the situation (see e.g. [3]). Notice that the predicted differential energy 
spectrum is generally flatter than the mirror matter case (or mini-charged dark matter 
case). Qualitatively, this can be very easily understood from the Eji dependence of 
the cross section, Eq.(ll). [The reason for this different Er dependence in the cross 
section is because of the fundamental interaction: WIMPs are weakly interacting and 
mirror matter is electromagnetic]. Evidently mirror matter works for two quite differ- 
ent reasons a) the cross section is electromagnetic which gives the differential energy 
spectrum the right shape and b) the value of fo(O') is much lower than the anticipated 

^Of course standard (SI) WIMPs are also disfavoured by several independent dark matter 
experinients[24, 25, 26] - quite unlike the case of mirror matter-type dark matterfl, 20]. 
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220 km/s due to the effects of the hght elements He' and/or H' (which dominates the 
mass of the halo) and also the mirror electrons (which contribute to the pressure) - 
this gives the spectrum the right normalization. 

In conclusion, we have further explored the mirror matter dark matter interpreta- 
tion of the DAMA/Nal experiment. Any proposed explanation of the annual modu- 
lation signal must also be consistent with the measured differential energy spectrum. 
It turns out that this requirement is quite non-trivial, nevertheless the simplest mirror 
matter interpretation of the DAMA/Nal annual modulation signal - a He' and/or H' 
dominated halo with small O' subcomponent - passes this test. It is striking that the 
simplest dark matter candidate coming from fundamental physics has just the right 
properties to fully explain the DAMA/Nal data. These are strong reasons to beheve 
that the dark matter problem has finally been solved. 
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Appendix 

Mirror matter-type dark matter galactic halo distribution 

Assuming a (standard) spherical isothermal dark matter halo, the condition of 
hydrostatic equilibrium gives: 

where P is the pressure (for a dilute gas, the pressure is related to the number density, 
n, via P = nT), p is the mass density and g{r) is the local acceleration at a radius r. 
This acceleration can be simply expressed in terms of the energy density, via: 

y = 1 r Gp47rr''^dr' . (13) 
Jo 

As is well known, the solution of Eq.(12,13) is n oc 1/r^, which gives the required flat 
rotation curve. 

The local velocity, Vrot{Ro) — 220 km/s at our location in our galaxy is related to 
the mass density via: 

v'rotm = ^ P^^r'dr. (14) 
Kq Jo 

If the mass of the galaxy is dominated by the halo component, then p oc 1/r^. The 
proportionality coefficient can be obtained from Eq.(14), which gives 

7;2 1 
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The gravitational acceleration at a radius, r, is then 

9 = \vlt ■ (16) 

Using this equation, we can solve the hydrostatic equilbrium condition, Eq.(12). As- 
suming the simplest case of a halo made up of a single particle species of mass m, 
Eq.(12) yields: 

T= — ^. 17 

2 ^ ' 

Defining the fo parameter such that the Maxwellian distribution function \j{y) — 
exp (^—^mv'^/T^] has the form 

f{v) = exp (^-^j (18) 

implies that Vq — 2T/m. Comparing this with Eq.(17) we find that in the simplest 
case of a single component halo, vq = Vrot (the standard result). 

Repeating this exercise for a spherical isothermal halo composed of several com- 
ponents, H',He',0'... at a common temperature T, we find pA' — Ma'Tia' oc 1/r^, 
with 



p^MpV, 



2 

rot 



(19) 

and 

vli^) ^ .20) 
v^t Ma'' 

In the above equation, fiMp is the mean mass of the particles comprising the mir- 
ror (gas) component of the halo {Mp is the proton mass). In determining the mean 
mass of the particles in the halo we must include the mirror electron component: at 
temperatures of order a keV, H', He' and other hght elements are fully ionized. 

Let us now consider two illustrative cases: a) the mass of the gas component of 
the halo is dominated by He' and b) the mass of the gas component of the halo is 
dominated by H'. In the first case, we have rie' = 2nHe' (since He' should be fully 
ionized), which means that iJ.Mp ~ Mhc'/^ and from Eq.(20) we have: 

v'^(He') 1 

-^—^ ~ -, [Case a) Mass of the halo dominated by He'] (21) 

In the second case, we have rie' = uh', which means that /xMp ~ Mh//2 and from 
Eq.(20): 

vi{H') 1 



,2 



[Case b) Mass of the halo dominated by H'] (22) 



'^rot 2 

In each case, the value of vo for a small O' component can be obtained from 



vo{0') - Vo{X')^Mx'/Mo' . (23) 

where X' = He' (first case) and X' = H' (second case) . [Of course the above equation 
is valid for any X'; it assumes only thermal equilibrium]. 
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